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Abstract 

Molten metallic nanoparticles have recently been used to construct graphene nanostructures with 
crystallographic edges. The mechanism by which this happens, however, remains unclear. Here, we 
present a simple model that explains how a droplet can etch graphene. Two factors possibly contribute 
to this process: a difference between the equilibrium wettability of graphene and the substrate that 
supports it, or the large surface energy associated with the graphene edge. We calculate the etching 
velocities due to either of these factors and make testable predictions for evaluating the significance of 
each in graphene etching. This model is general and can be applied to other materials systems as well. 
As an example, we show how our model can be used to extend a current theory of droplet motion on 
binary semiconductor surfaces. 

PACS numbers: 81.16.Rf, 68.65.Pq, 61.48. Gh 



1 Introduction 

Graphene-derived nanomaterials are emerging as a unique class of low-dimensional structures, in part 
due to their remarkable electronic properties and intriguing chemical reactivity [TJ [2]. Graphene's two- 
dimensionality makes it compatible with current planar device architectures. In addition, its electronic 
properties can be a strong function of how it is structured [3l |U [5j [6l [7]. Because of these characteris- 
tics, a significant amount of research is currently focused on constructing graphene-based nanoelectronic 
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devices, which hold enormous promise for high-performance post-silicon applications. Desirable devices in- 
clude nanoribbons with crystallographic edges, which would potentially combine the exceptional electronic 
properties of graphene with rational control of device characteristics by modification of the graphene archi- 
tecture. 

Experimentally, constructing these devices is a formidable challenge [2J. Most current nanolithography 
techniques produce graphene structures with rough, noncrystalline edges, which can be quite detrimental to 
device performance [3J [31 [5] GO EH HU HI]- Alternative methods are being explored, although a good deal 
of work is required to optimize and extend them to practical usage. In previous work, my colleagues and I 
developed a technique by which molten metallic nanoparticles etch straight trenches in few-layer graphene 
flakes supported on amorphous SiC>2 |13) . Strikingly, we found that these trenches span lengths of > 1/zm 
along directions commensurate with the graphene honeycomb lattice and extend down to the underlying 
substrate through the graphene. Further work has shown the edges of these trenches to be crystallographic 

The microscopic mechanism by which this crystallographic etching occurs is unclear. Here, we propose 
a simple model for this process, motivated by the results of etching experiments as well as theoretical 
models for "running" droplets [151 116] . Two factors potentially drive crystallographic etching by moving 
droplets: a difference between the equilibrium wettability of the graphene and the underlying substrate, 
or the high surface energy associated with a reactive graphene edge. The analysis presented here yields a 
number of testable predictions that could be used to characterize the significance of each of these factors 
in controlling graphene etching. Furthermore, the ideas presented here are relevant to a number of other 
materials systems as well. As an example, we study the case of droplets formed during the evaporation of 
binary semiconductors. We use our model to show how a current theory of this process may be extended to 
better agree with experimental observations. 

2 Model and Results 
2.1 Experimental Motivation 

As reported previously [13] , a commonly- used preparation process for etching of few-layer graphene samples 
supported on Si02 involves uniformly coating them with a solution of iron salt. These samples are then 
heated in a hydrogen/argon environment at 900°C for a given period of time, typically on the order of 1 
hour. Under these conditions, the iron salt forms nanoparticles with diameter ~ 20nm [17] . Because this 
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temperature is above the melting point of the nanoparticles, they form molten droplets [T51QI5]. Driven by 
Brownian motion, these diffuse over the surface of the amorphous SiC>2 substrate or the graphene before 
encountering a graphene edge. 

Because they are not part of the bulk sp 2 -bonded network, carbon atoms at graphene edges are chemically 
reactive [201 HTJ [22] . (On the other hand, atoms not at edges are inert and likely do not react with the molten 
droplets. This is supported by observations that etched trenches begin only at graphene edges.) These atoms 
likely undergo catalytic hydrogenation, a chemical reaction catalyzed by the molten Fe droplet in which the 
carbon atoms react with the hydrogen supplied from the surrounding environment, forming methane gas 
[13) . This process is thus initiated when the thermally-driven droplet randomly encounters such an edge. 
The droplet removes the contacting atoms from the graphene sheet, forming a trench, and propagates the 
reactive edge further into the graphene sheet - we term this the step edge. This is shown schematically in 
Figure 1. Experimental observations show these trenches to be straight over micrometers, thus indicating 
that lateral motion of the droplet due to thermal effects can be neglected while considering the etching 
process. Thermal effects may become significant for droplets much smaller than those considered here. 

This procedure has since been extended to the etching of monolayer graphene, the surface of bulk highly- 
oriented graphite, and carbon nanotubes using a variety of metal catalyst materials and at a range of high 
temperatures [231 H2 (HI (23 [2H1 [23 12H1 122] • The temperatures used in these experiments are similar to or 
higher than the nanoparticle melting temperature, which is well known to be significantly lower than the 
melting temperature of the bulk material [30| . Thus, the nanoparticles are likely to either be fully molten or 
have a molten surface layer. Furthermore, the etched trenches and the manner in which they form appear to 
be similar in these experiments. Taken together, these observations indicate that metal nanoparticle etching 
of graphene at high temperatures may have a common origin, and a model treating the nanoparticles as 
chemically-reactive droplets may capture the underlying physics of this process. Two criteria guide the 
continued etching along straight crystallographic lines: the choice of the initial direction of the droplet and 
the constant velocity of the droplet along that direction. 

2.2 Outline Of The Model 

The choice of the initial direction of the droplet as it meets the graphene edge is a result of the chemistry of the 
etching process. Experimental observations indicate that the etched trenches extend down to the underlying 
substrate |13] , suggesting that removal of carbon is fast compared to the motion of the droplet. Experiments 
on bulk samples [31] as well as density-functional theory simulations [14] indicate that this process is likely 
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to be strongly anisotropic. That is, the energy required for removal of carbon atoms from a graphene sheet 
by catalytic hydrogenation is significantly reduced along particular crystallographic directions. Thus, when 
a molten droplet meets a reactive graphene edge, etching occurs fastest along whichever crystallographic 
direction requires the least amount of energy for carbon removal to occur. While this step is not the focus 
of this paper, we stress that a deeper understanding of it is crucial to a full understanding of nanoparticlc 
etching of graphene. 

Once the carbon atoms in contact with the molten droplet are removed and the initial droplet direction 
is chosen, the droplet moves along the crystallographic direction. The droplet must continue along that 
direction for successful etching to occur - this step is the main focus of this paper. We model the molten 
nanoparticle as a liquid droplet of characteristic size r, density p, and liquid- vapor surface tension ji v . For 
the iron nanoparticles used in graphene etching, r « 20nm 13 . Published values of the density and surface 
tension give p w 7g/cm 3 and ji v « 2N/m [33J [33]. These values yield a capillary length A = \/ji v /pg ~ 
5mm, where g is gravitational acceleration - because r/X <g; 1, we can safely assume that gravitational effects 
are negligible. Experimental etching velocities are on the order of nanometers per second [13j . implying that 
the particle Reynolds number Re = poVr/r/o <C 1, where po and 770 & re the density and viscosity of air and 
V is the etching velocity. The minimal model presented here focuses on key details of graphene etching 
by nanoparticles to obtain a qualitative understanding of this process. Specifically, we derive the existence 
of a force driving droplet motion in one direction and investigate the dependence of the droplet velocity 
on relevant physical parameters. Potential extensions of this model could explicitly treat hydrodynamic 
flows within the droplet, the atomistic details of the graphene edge, and the effects of non-idealities such 
as contact-line pinning and hysteretic effects resulting from surface heterogeneities. Furthermore, one could 
incorporate electrostatic effects that may influence droplet characteristics [34] . 

At a given time after the choice of the initial etching direction, the droplet straddles the step edge. For 
simplicity, we consider a two-dimensional cross-section, as shown in Figure IB. The advancing side is on 
the graphene at point g with advancing contact angle 6g, while the receding side of the droplet lies on the 
substrate at point s with receding contact angle Q r s . The equilibrium contact angles of the droplet deposited 
fully on graphene or on the substrate are defined to be 9 g and 6 S , respectively. The graphene- liquid, graphene- 
vapor, substrate-liquid, and substrate- vapor interfacial energies are 7^ , j gv , 7^ , and j sv , respectively. The 
thermodynamic driving force for droplet motion is a result of the unbalanced Young's force at the contact 
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line, and so the total driving force on the droplet goes as 

Fe ~ r • [(<y gv - 7 9i ) - (j sv - j si )} (1) 
= r- 7to (cos^-cos^) (2) 

There are thus two factors that will drive a droplet to keep moving in one direction: (1) a difference between 
the equilibrium wettability of graphene and of the substrate that supports it, in which case 9 g and 9 r a are 
determined primarily by 9 g and 9 Sl or (2) the large surface energy associated with the graphene step edge, in 
which case 9 a g and 9 r s are determined primarily by the magnitude of the energy cost associated with having 
a step edge. We first focus on case (1) and ignore the structure of the graphene step edge. 

2.3 Motion Driven By Differing Wettabilities 

As the droplet straddles the graphene step edge, it will remove the underlying reactive carbon via catalytic 
etching. To relate the thermodynamic driving force to the nature of the underlying substrate, we adapt a 
model proposed by Brochard and de Gennes [15j . For the given droplet configuration, we assume an areal 
density of carbon atoms that monotonically varies from 1 at the droplet advancing side to a value <p < 1 
at the droplet receding side. In general, we expect the surface energies J s i,"fsvi"fgii and "f gv to depend on 
4>. Specifically, we expect the difference between (7^ — j sv ) and (7^ — J gv ) to increase with decreasing 
4>. To first order in (j>, this is manifested in the relation (7^ — 7 S „) — ("f g i — J gv ) = 7o(l — <fi), where 70 
is a parameter that reflects the increase in the differential surface energy for a reaction that has achieved 
completion. Relating this to the equilibrium contact angles via the Young equations j gv — "f g i = 7;„cos 9 g 
and 7 S1) — 7^ = 7;„cos 9 Sl we obtain 

ji v (cos 9 g - cos 9 S ) = 70(1 - 4>) (3) 

Assuming that the equilibrium wettabilities of the graphene and of the substrate are the prime determinant 
of 6g and 6 r s , as in case (1), we consider the limit where 9 g « 9 g and 9 r s « 9 S . As a result, Fg/r w 70(1 — cf>). 
The chemistry of the etching process determines (f>. Assuming first-order reaction kinetics, (f> — e~ r / VT , 
where V is the velocity of the droplet and r is the characteristic time constant of the reaction. This yields 
a direct relationship between the dynamics of the etching reaction, quantified by the time constant t, and 
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the thermodynamic driving force for droplet motion, Fg: 

F e ~r- 7o (l - e^'^) (4) 

We treat the time-averaged droplet motion as a damped response to the force driving droplet motion, 
balancing Fg with a drag force F^ — £V, where £ is the drag coefficient. Equating these forces results in a 
transcendental equation for the droplet velocity: 

V ~ a{l-e~ p/v ); a = 7o r/£, /? = r/r (5) 

Crucially, the droplet velocity depends on two characteristic velocities, a and /3, or what we term the driving 
velocity and the reaction velocity, respectively. The numerical solution to Equation 4 is presented in Figure 
2. The droplet velocity increases monotonically with increasing a or increasing f3. 

We make two assumptions to directly relate a to experimentally measurable quantities. First, 9 g — > 9 g 
and 9 r s — > 9 S in the limiting case of <j> ~ Equation 2 then gives 70 = 7i!,(cos 9 g — cos 6 S ). Second, 
we explicitly consider the nature of the drag coefficient £ in Equation 4. This can generally be written as 
£ = Wf( r ■>&*)■) where the dynamic contact angle 9* = cos -1 [(cos 9 g + cos 0£)/2], / is a function of r and 
9* appropriate to the system under consideration, and 77 is the droplet viscosity. Because we assume in this 
section that equilibrium wettabilities are the strongest determinant of the advancing and receding contact 
angles, we take 9* « cos -1 [(cos 9 g + cos 9 s )/2}. f(r,8*) now remains to be determined. 

Following previous work |35j . we assume that dissipation occurs primarily via viscous processes at the 
droplet-substrate contact line, as is typically assumed when considering low Reynolds number free surface 
flows with a moving contact angle. In particular, we apply the lubrication approximation to a "wedge" of 
angle 9*, assuming a two-dimensional Poiseuille velocity field in the wedge [55]. This yields f(r,9*) oc 1/9*, 
where the geometrical details of the system are lumped into I = hi\x max /x m i n \, & logarithmic function 
of two parameters that describe the geometry of the region within which dissipation is concentrated. 
Xmax is typically on the order of the droplet size (« 10 — lOOnm), while x m i n is typically taken to be 
on the scale of individual molecules (« 0.1 — lnm), giving I ~ 1 — 10. Combining these results gives 
£ oc 7/r/cos -1 [(cos 9 g + cos 9 s )/2]. 

We need to consider the validity of the lubrication approximation, which strictly speaking only applies to 
small values of 9* , in the situation presented here for a broad range of contact angles < 180°. For example, 
dissipation may cease to be localized at the droplet-substrate contact line for sufficiently large values of the 
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contact angle. However, because of the large viscosity of the droplet relative to the medium surrounding 
it, it is unlikely that the scaling behavior will change in the situation considered here. Indeed, the more 
general case of larger 9* values has been considered theoretically. While the ~ 1 /8* scaling of the dissipation 
presented here crosses over to a ~ 9* 2 scaling with increasing contact angle, it is important to note that 
this crossover occurs when 9* = (?7o/??) _1/ ' 3 [SZ]- Using Sutherland's empirical formula [35] to estimate the 
viscosity of the air surrounding the droplet gives ryo ~ 0.02 — 0.06mPa-s, while literature values of liquid 
metal viscosities are in the range 77 ~ 2 — 6mPa-s. These parameters yield a crossover contact angle > 180°, 
indicating that the ~ 1/9* scaling of dissipation considered here is appropriate for our system. We thus have 
the characteristic velocities explicitly in terms of the experimentally-measurable equilibrium contact angles: 

a cc 7;„?7 -1 (cos 9 g — cos 9 S ) ■ cos^ 1 [(cos 9 g + cos 9 s )/2] , j3 = r/r (6) 

These results have important implications. As shown in Figure 3, the sign and magnitude of the droplet 
velocity is strongly dependent on the values of 9 g and 9 S . Notably, V > only when 9 g < 6 S . While only a 
few clear equilibrium contact angle measurements of 9 g and 9 S exist for typical nanoparticle materials, these 
data suggest that the mechanism driving droplet motion outlined in case (1) may play a significant role in 
graphene etching, as shown in Figure 3. More details are provided in Appendix A. 

2.4 Motion Driven By Edge Energetics 

The structure of the graphene step edge may also contribute to driving droplet motion in graphene etching. 
This driving force is the dominant mechanism of case (2), in which the equilibrium wettabilities of graphene 
or the substrate do not play a significant role. Instead, the difference between the surface energies j gv and 
7 st , (and hence, the values of 9 g and 9 r s ) will be determined primarily by the energetics of the step edge. We 
will consider this case using a treatment similar to that used in Section 2.3. 

Because the atoms at a freshly-created step edge have dangling bonds, we associate a surface energy with 
each step SOI HD S3- We make the simplifying assumption that j g i — 7 S ; and ^ gv = A + 7 S „, where A 
is the total surface energy density associated with the overall step structure across the droplet profile. This 
step structure is shown schematically in Figure 1. Thus, as given by Equation 1, the driving force for droplet 
motion goes as Fg ~ r • A. We adopt the previous discussion to this situation: coarse- graining over the 
atomic details of the step edge, we assume an areal density of steps that monotonically varies from 1 at the 
droplet advancing side to a value <j> < 1 at the droplet receding side, such that A = Aq(1 — <t>). The results of 



7 



the previous case thus carry over, with the key replacement of the phenomenological surface energy 70 with 
the step energy Ao- The droplet velocity again depends on the two characteristic velocities a and j3 as given 
by Equation 4. As in case (1), the droplet velocity increases monotonically with increasing a or increasing 
ft, as shown in Figure 2. However, a is now determined by different physical quantities as manifested in the 
new relation a = Aor/£. 

We can relate the step energy to the structure of the step [3H1 HOI Ell HH S3] ■ In particular, assuming that 
the energy of an individual step is proportional to the number of atoms making up the step, Ao oc (t/a)-Eb on d 
where t is the step height, a is the average bond length of the graphene in the direction normal to its surface, 
and Ebond is the energy required to form one dangling bond, with the proportionality constant incorporating 
the geometric details of the step structure. 

Unlike the previous case, we cannot express the advancing and receding contact angles 8® and 9 r s , and 
thus the drag coefficient £, in terms of experimentally-determinable equilibrium contact angles. However, as 
in the previous case, we can assume that dissipation occurs primarily at the droplet-substrate contact line. 
This gives £ oc ?yr, and hence 

a oc (t/a) ■ E bond r]" 1 , fi = r/r (7) 

Finally, we note that because the catalytic hydrogenation reaction occurs very quickly, as noted in Section 
2.2, the reaction time constant r is likely to be small relative to the time the droplet takes to move its own 
length. It is thus probable that the etching process is not limited by the speed of the reaction and a is 
the key determinant of the droplet velocity for both of the driving forces described here. Furthermore, this 
supports our assumption that the metal catalyst is not saturated with an appreciable amount of carbon 
during the etching process and does not lose catalytic behavior as the reaction progresses. 

3 Discussion 

3.1 Relevance To Graphene Etching 

We have focused our discussion on developing predictions for the time-averaged etching velocity, V. Not 
only is maximizing V crucial to optimizing this etching process, but it is also likely to be the most easily 
measurable physical parameter in experiments. Measuring the etching velocity can be done in a number of 
ways. One possibility is to use real-time transmission electron microscopy [44] . but this can be experimentally 
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challenging. Another approach is to pattern the catalyst particles at a graphene step edge [57]. Measuring 
the length of the trenches formed after a given etching time directly yields the distribution of time-averaged 
etching velocities. 

The force driving droplet motion in graphene etching is likely to have contributions from both factors 
described in Sections 2.3 and 2.4 - namely, (1) a difference between the equilibrium wettability of graphene 
and of the substrate that supports it and (2) the high surface energy associated with forming a graphene 
step edge. An outstanding question is whether one of these is a stronger determinant of droplet motion 
than the other. Crucially, different physical parameters determine the contributions from factors (1) and 
(2), enabling their importance to be tested experimentally. 

If factor (1) plays a significant role, etching velocities will be a strong function of the wettability of the 
underlying substrate, as indicated in Equation 5. This effect can be tuned by studying etching of the same 
material on entirely different substrates. To assess the experimental feasibility of this test of factor (1), we 
list in Table 1 literature values of 9 S for Fe, Ni, Co, and Ag droplets on a number of materials that are 
commonly used as substrates for graphene |451 1461 147) . Comparing these to the data presented in Figure 3 
suggests that if factor (I) influences the motion of droplets, the etching velocities for Fe should follow the 
trend vsio 2 > v ai 2 o 3 > vsi 3 Ni > i>sic, where Vi indicates the etching velocity on substrate i. Furthermore, 
the etching velocities of Ni or Co on SiC>2, AI2O3 and SisN4 should be similar to each other, and significantly 
greater than the etching velocities of Ni or Co on SiC. The case of Ag is a unique one. Because we do not 
expect factor (1) to contribute to its etching behavior, the etching velocities of Ag should be independent 
of the substrate used. A more practical alternative is to study etching on a substrate that has been surface- 
treated to various degrees, for example, by coating it with different concentrations of a suitable chemical 
species. The etching velocities will be a strong function of the wettability of the graphene itself, as well. The 
graphene wettability can be varied via suitable surface treatment, such as by using a range of strengths of 
oxygen plasma etching [48] or by sonicating the graphene in different solvents [49] . 

If factor (2) plays a significant role, etching velocities will be a strong function of the thickness of the 
graphene. A systematic study of etching velocity as a function of graphene thickness remains to be done. 
The thickness is easily varied experimentally. For example, micromechanical cleavage of highly-oriented 
pyrolytic graphite using Scotch tape can result in monolayer graphene flakes and few-layer graphene flakes 
with thicknesses up to tens or hundreds of layers. As described in Section 2.4, the etching velocity should 
increase linearly with the thickness of the graphene. Performing this study using Ag droplets would be 
particularly interesting. Because we have ruled out the possible contribution of factor (1) to etching by Ag, 
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as described in Appendix A and Figure 3, we expect factor (2) to be the only driving force for etching. 



3.2 Relevance To Binary Semiconductors 

"Running" droplets similar to those considered here form on the surface of evaporating binary semiconduc- 
tors. Unlike the inert substrate in the case of graphene etching, the underlying semiconductor substrate 
acts as a reservoir for droplet formation. Experimental advances have led to the direct imaging of droplet 
formation and movement. As recently shown by Tersoff et al. [50] . droplets form on the surface of the 
binary compound when it is annealed at high enough temperature, move in directions commensurate with 
the underlying crystal lattice, and leave behind a trail potentially similar to the etched trenches described 
earlier for graphene etching 51 . To explain this, they model this process by considering the thermodynamics 
of noncongruent evaporation, in which one chemical component of the binary compound evaporates slower 
than the other. We believe that this model is incomplete and propose that the ideas presented in Section 
2.4 can provide a more complete physical picture. Here we briefly summarize Tersoff et aVs argument and 
extend it using our model. 

For temperatures T higher than T c , the upper limit for congruent evaporation, Ga forms droplets on the 
GaAs crystal surface because the Ga chemical potential of the surface /j, s is larger than the liquidus chemical 
potential fii. If one assumes the GaAs crystal surface to be totally isotropic, the driving force for droplet 
motion presented in Equation 1 is zero, since ^ gv = 7 S1 , and 7 g ; = 7^. Tersoff et a/.'s model breaks this 
symmetry by assuming that while the GaAs surface is topographically unchanged, a droplet displacement to 
one side (e.g. by thermal fluctuations) will cause the surface energy of the newly-exposed part of the GaAs 
surface to be at 7 OT = 7 s „(n . On the other hand, the newly-covered part of the GaAs surface will still remain 
at 7 gl , = 7 S „( S ) ^= Jsvtl)- Because the surface and liquid both contain Ga, 7 g ; is assumed to be equal to 7^. 
Equation 1 then yields Fg ~ r ■ (p/ gv — j sv ) = r ■ (7 S „( S ) — J sv a))- Tersoff et al. relate this to the difference 
between the Ga chemical potential of the newly-covered surface fi s , which is a function of T, and the Ga 
chemical potential of the newly-exposed surface which is taken to be the liquidus chemical potential. In 
particular, they use an equilibrium expansion of the chemical potential around /ii to find 



where E is the Ga surface formation energy and n is the excess Ga surface density. They then assume 



Fg ~ r ■ (7 sl) ( s ) - 7at)(i)) 



(8) 




(9) 
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Ha — M oc T — T c to derive their main result, namely that the driving force for droplet motion varies 
quadratically with temperature around T c . The strength of this driving force is determined by the pre- factor 
in Equation 7. Unfortunately, the value of this pre-factor cannot be explained by any current theory. Indeed, 
there appears to be very little discussion in the literature investigating the variation of E and n with fi 52 . 

While this model sets important groundwork for understanding running droplets on binary semiconduc- 
tors, we believe it is incomplete. In particular, it neglects two crucial experimental observations. First, the 
model does not incorporate any directionality in the droplet motion. Because the droplet motion is assumed 
to be started by thermal fluctuations, the direction of the droplet trajectory should be randomly oriented. 
This prediction contradicts experimental observations, in which droplets are seen to clearly move preferen- 
tially along crystallographic directions [50] . Second, the model does not incorporate topographic changes in 
the semiconductor surface. Specifically, observations of droplets moving on GaAs |50j and on other binary 
compounds [51] have shown that trenches or step structures possibly form in the surface as droplets move 
over it. The steps are typically oriented normal to the direction of droplet motion |51j . similar to the step 
structure discussed in Section 2.4 for graphene and schematically shown in Figure 1. 

Taken together, we believe these observations as well as our model for graphene etching imply a possible 
modification to Tersoff et aVs model. This extended model may better capture all the experimental obser- 
vations for running droplets on binary semiconductor surfaces. We propose that surface reconstruction, such 
as the formation of trenches and steps in the GaAs surface, occurs fastest along specific crystallographic 
directions, as is common for many crystals. This process sets the initial direction of droplet motion - not 
thermal fluctuations. Furthermore, we use the model of Section 2.4 here to determine the droplet velocity. 
This gives V ~ a(l — e~^l y ) with a = Aor/£ and j3 = r/r, where r is related to the rate at which Ga is 
removed from a step. Unlike the graphene etching case, because the GaAs substrate acts as a reservoir for 
the Ga droplet, Ao is not necessarily associated with the density of dangling bonds associated with a step 
edge as assumed in Section 2.4. Rather, we hypothesize that this phenomenological parameter incorporates 
the temperature dependence of the Ga droplet velocity. For example, we suggest that Ao oc (/x s — /^z) 2 , in 
agreement with physical intuition: when /i s — fii, Ga atoms at the GaAs step edge will not have an energetic 
preference to be in either the droplet or in the GaAs surface. This represents an effective GaAs edge energy 
Ao that is zero at /i s , and increases for variations of fj, s around /i/ . Motivated by the discussion of Tersoff et 
al., we assume that fi s — Hi oc T — T c . The treatment presented here thus yields a droplet velocity V that 
increases monotonically with (T — T c ) 2 , consistent with the data presented in [50] . 
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4 Conclusion 



The development of a straightforward technique to carve out straight, large-scale features with crystal- 
lographic edges in graphene would be invaluable. For example, this would represent a way of achieving 
graphene devices with tunable electronic band structures and exceptional performance. Recent experiments 
using molten metallic nanoparticles to construct graphene nanostructures with crystallographic edges arc 
one step toward this goal, and as such, understanding the mechanism by which this happens is crucial. The 
simple model presented here gives a qualitative understanding of two factors driving droplet motion and 
etching: (1) a difference between the equilibrium wettability of graphene and the substrate that supports 
it and (2) the large surface energy associated with the graphene step edge. It is unclear whether both of 
these factors play a role in driving etching. Crucially, the time-averaged etching velocities resulting from 
each of these factors depend on different physical parameters. The model presented in this paper makes 
clear predictions for each of the factors that can be tested, and we anticipate that the results presented here 
will help guide experimental efforts using nanoparticles to etch graphene. Furthermore, we propose that this 
model is applicable to other materials systems. As an example, we apply it to the case of running droplets 
formed during the evaporation of binary semiconductors. We use our model to extend a current theory of 
this process, showing how the work presented in this paper can be used to improve the agreement of this 
theory with experimental observations. 
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Table 1: Literature values of 9 g and 9 S for droplets of different materials on a variety of substrates at high 
temperatures. 
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A Comparison To Previous Experiments. 

To assess the applicability of the results presented in Section 2.3 - in particular, Equation 5 - to graphene 
etching, we examined the literature for contact angle measurements for nanoparticle materials relevant to 
experiments. These are Fe, Ni, Co, or Ag. The results are given in Table 1. Systematic contact angle 
measurements of molten nanoparticles on single- or few-layer graphene have, to our knowledge, not yet 
been performed. Thus, to get a sense of likely values of 9 g , we turn to high-temperature measurements 
on bulk graphite substrates. Indeed, recent work suggests that the thickness of the graphene may not 
play a significant role in altering its wettability |48j . These data do not appear to be strongly sensitive to 
temperature, so we expect them to be similar to the actual contact values in graphene etching. 

The data presented in Table 1 are typically for "clean" droplets on the respective substrates. Various 
chemical processes might change these values. For example, significant concentrations of carbon saturated 
in the metal droplet are known to modify the value of 9 g . However, values for "clean" unsaturated droplets 
are probably most relevant for etching. Removal of carbon atoms via catalytic hydrogenation occurs quickly, 
as evidenced by observations that etched trenches extend down to the underlying substrate. Thus it is quite 
likely that the metal catalyst at a step edge is not saturated with an appreciable amount of carbon during 
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the etching process. 

As discussed in Section 2.3, V > only when 6 g < 9 S . The data presented in Table 1 suggest that 
graphene etching by Fe, Ni or Co satisfy this condition, and thus the mechanism discussed in Section 2.3 
may play a role in driving the motion of these droplets. On the other hand, graphene etching by Ag does not 
satisfy this condition. Crucially, this implies that etching by Ag is driven solely by the mechanism presented 
in Section 2.4 or by some other mechanism not considered in this paper. 
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Figure 1: Schematic of the experiment. (A) Top-down view of graphenc etching by a molten nanoparticle, 
showing possible step structure of graphene layers at the edge undergoing etching. (B) Two-dimensional 
representation of etching, showing advancing and receding contact angles. The grey underlying layer is the 
substrate while the black line being etched away is the graphene, showing an example of a step edge. These 
diagrams are not to scale. 
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Figure 2: Droplet velocity. Dependence of the droplet etching velocity V as a function of the characteristic 
velocities a and /?, as given in Equation 4 of the text. 
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Figure 3: Influence of substrate and graphene wettability. Wettability "state diagram", showing 
droplet etching velocity V as given in Equation 5 of the text. Shaded region represents positive values of 
V, as given by the color scale and contours. Rectangular regions and unfilled circles indicate where different 
nanoparticle materials lie on this diagram. 
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